Hadronic molecule structure of the Y"(3940) and F(4140) 
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We report on further evidence that the Y(3940) and the recently observed Y(4140) are heavy 
hadron molecule states with quantum numbers J PC = ++ . The Y(3940) state is considered to 
be a superposition of D* + D*~ and D*°IF°, while the Y(4140) is a bound state of D* + and D* 
mesons. For the first time we give predictions for the strong Y(3940) — > J/ipiv, Y(4140) — > J/ip4> 
and radiative T(3940)/T(4140) — > 77 decay widths in a phenomenological Lagrangian approach. 
Results for the strong J/ipV (V = u, <f>) decays clearly support the molecular interpretation of the 
T(3940) and Y(4140). The alternative assignment of J PC = 2 ++ is also tested, giving similar results 
for the strong decay widths. 

PACS numbers: 12.38.Lg, 12.39.Fe, 13.25.Jx, 14.40.Gx 
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The recent announcement [l[ of the narrow state 
Y(4140) by the CDF Collaboration at Fermilab raises 
the prospects of possibly uniquely identifying the struc- 
ture of a meson resonance which does not fit in the 
conventional quark-antiquark picture. This latest re- 
port of a narrow structure with charmonium-like de- 
cay modes is a continuation of previous discoveries of 
such states which are not easily explained as quark- 
antiquark configurations. Possible alternative interpre- 
tations involve structures such as hadronic molecules, 
tetraquark states or even hybrid configurations (for re- 
cent reviews see e.g. Refs. [lli]). Now the CDF Collab- 
oration has evidence of a narrow near-threshold struc- 
ture, termed the Y(4140) meson, in the J/tp(f> mass spec- 
trum in exclusive B + — > J /ip<j)K + decays with the mass 
m Y {iiw) = 4143.0 ± 2. 9(stat) ± 1.2(syst) MeV and nat- 
ural width r y(4140) = 11.7i|;g(stat) ±3.7(syst) MeV Q. 
As already stressed in the new structure Y(4140), 
which decays to J/ip<fi just above the J/ip<f> threshold, 
is similar to the previously discovered Y(3940) 0, @], 
which decays to J/ij)ui near this respective threshold. 
The mass and width of the F(3940) resonance are: 
TOy(3940) = 39 43 ± ll(stat) ± 13(syst) MeV, Tyamo) = 
87±22(stat)±26(syst) MeV (Belle Collaboration f|) and 
my (39 4o) = 3914.6li;|(stat) ± 2.0(syst) MeV, r y(3940) = 
34±f (stat) ± 5(syst) MeV (BABAR @|). Both observed 
states, V(4140) and Y(3940), are well above the thresh- 
old for open charm decays. A conventional cc charmo- 
nium interpretation is disfavored, since open charm de- 
cay modes would dominate, while the J/i>4> or J/ipu de- 



cay rates are essentially negligible 0, 0]. Note, current 
data imply a lower bound of r(Y(3940) — > J/ipu) > 
1 MeV [4j, which is an order of magnitude higher than 
typical rates between known charmonium states. This 
could be a signal for nonconventional structure of the 
Y(3940). As a first follow-up to the CDF result, it is 
suggested in Q that both the V(3940) and V(4140) are 
hadronic molecules (i.e. bound states of mesons induced 
by the strong interaction). These hadron bound states 
can have quantum numbers J PC = ++ or 2 ++ whose 
constituents are the vector charm D*(D*) mesons: 



|Y(3940)) 



V2 

|V(4140)) = \D*+D i 



(\D*+D*-) + \D*°D* )) , 



(1) 
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The authors of Ref. (jj show that binding of the above- 
mentioned meson configurations can be achieved in the 
context of meson-exchange potentials generated by the 
Lagrangian of heavy hadron chiral perturbation theory 
(HHChPT) 0-[I(l. Earlier results based on the pion- 
exchange mechanism already indicated that the D*D* 
system can form a bound state ll| . Binding in the D*D* 
channel can be induced by 77 and 4> meson exchange 0]. 
A first QCD sum rule study cannot support the claim 
that the D*D* system binds [12}. This issue remains to 
be studied. 

In this Letter we report on a first quantitative predic- 
tion for the decay rates of the observed modes V(3940) — > 
J/ipu> and V(4140) — * J/tp(f> assuming the hadronic 
molecule structures of Eq. |T|) with quantum numbers 
J PC = ++ . Results will be shown to be fully consis- 
tent with present experimental observations, strengthen- 
ing the unusual hadronic molecule interpretation. Fur- 
ther predictions are given for the radiative two-photon 



2 



decays of these states. Finally, we also consider the al- 
ternative J PC = 2 ++ assignment for the Y states. 

The method of determining the decay rates is based 
on an effective Lagrangian which includes both the cou- 
pling of the molecular-bound state to their hadronic con- 
stituents and the coupling of the constituents to other 
hadrons and photons. In Refs. [l3| we developed the for- 
malism for the structural study of other recently observed 
meson states (like D* (2317), D s i(2460), V(3872), •••) 
as hadronic molecules. The composite (molecular) struc- 
ture of the V(3940) and F(4140) states is defined by 
the compositeness condition Z = (see also 

Refs. [l3j). This condition implies that the renormal- 
ization constant of the hadron wave function is set equal 
to zero or that the hadron exists as a bound state of its 
constituents. Decay processes are then described by the 
coupling of the final state particles via one-loop meson 
diagrams to the constituents of the molecular state (see 
details in [l3j]). 

For the observed V(3940) and F(4140) states we adopt 
the convention that the spin and parity quantum num- 
bers of both states are J PC = ++ . Presently, except for 
C = +, the J p quantum numbers are not unambiguously 
determined yet in experiment. For example, the V(3940) 
is also discussed as a J PC = 1 ++ charmonium candi- 
date [J], but ++ is not ruled out. Their masses are ex- 
pressed in terms of the binding energy ey as 777,y( 3 g 4 o) = 
2m D , - ey(394o) and m Y (4i40) = 2m D « - £y(4i40) > where 
771 £)» = m D ,+ = 2010.27 MeV and ttid* = m D *+ = 
2112.3 MeV are the masses of the constituent mesons. 
Since the observed masses are relatively far from the cor- 
responding thresholds we do not include isospin-breaking 
effects (i.e. we suppose that charged and neutral non- 
strange D* mesons have the same masses). Following 
Ref. [7J we consider the V(3940) meson as a superposi- 
tion of the molecular D* + D*~ and D*°D*° states, while 
the F(4140) is a bound state of D* + and D*~ mesons 
(see Eq. (JTJ). The coupling of the scalar molecular states 
to their constituents is expressed by the phenomenologi- 
cal Lagrangian: 



C Y (x) = g Y Y ij {x)Jy i Ax) 



(2) 



where gy is the coupling constant; Yy is the 3x3 matrix 
containing a nonet of possible hidden and open flavor Y 
states which can be composed of vector D*(D*) mesons: 







Y — 



Y„ 



Y, 



K* 



s/2 V2 p 

Y° Y 

Yp "71 Vl K 
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yO 
1 K- 



(3) 



In addition to the detected Y u = Y"(3940) and Y$ — 
y(4140) states, one can propose an isotriplet of non- 
strange states Y+ = (D*+D*°), Y~ = {D*-D*°), 
Y° = {D*°~D*v - D*+D*-)/y/2 and two isodoublets 



Y o 

K" 
)*+ i 



(D*-D*+) 
). We ex- 



of strange states Yi, = (D*°D 
and Y K , = (D*°D*-), Y°~ = (D*+D* 
pect that the masses of the , Y® states are close 
to the V(3940) mass, while the masses of the other 
four states Y^, 



YE, , Yjl, could be approximately 4040 



K , 

MeV ~ mo* + mr>* — 80 MeV (a typical value for 
the binding energy, as in the case for the V(3940) and 
V(4140) states). In analogy with the V(3940) and 
Y(4140) states, we suggest that the new hypothetical 



states Yp ±(0) (3940) and i£. w (4040) can decay into J/ifjp 
and J/tpK* pairs, respectively. Note, we use the nota- 
tion Yy for the nonet of Y states since it decays into 
J/tpV pairs. Jy ti of Eq. |(3J) is the current composed 
of the constituents of the respective hadronic molecule 
Yij. The simplest form of the hadronic currents J Yij 
is given by: Jy tj (x) = g^ v J d 4 y$(y 2 )J^ (x,y), where 



,±(0) 



j)1U 



(x,y) = D*/(x + 



%)D* u \x - § ). Here $(y 2 ) is the 
correlation function describing the distribution of the 
constituents inside the molecular states Y . For sim- 
plicity we adopt a universal equivalent function for all 
states. A basic requirement for the choice of an explicit 
form of the correlation function $ is that its Fourier 
transform vanishes at a sufficient rate in the ultravio- 
let region of Euclidean space to render the Feynman di- 
agrams ultraviolet finite. We use a Gaussian form of 
$(y 2 ): $(p|/A^) = exp(-p|/Ay), where p E is the 
Euclidean Jacobi momentum. Here, Ay is a size pa- 
rameter with a value of about 2 GeV - a typical scale 
for the masses of the constituents of the Y states. The 
coupling constants gy are determined by the composite- 
ness condition [l^]- 17|: Zy = 1 — Sy(rriy) = 0, where 
T, Y (m Y ) = dT l y{p 2 ) j dp 2 \ p 2 =m 2 r is the derivative of the 
mass operator Ey generated by Cy(x). 

To determine the strong Y — > J /ipV and two-photon 
Y — > 77 decays we have to include the couplings of 
D*{D* S ) mesons to vector mesons (J/tp, oj, (f>) and to 
photons. The couplings of J/tp, w, (j) to vector D*(D*) 
mesons are taken from the HHChPT Lagrangian Hoj: 



-D'D'V 



D 7 



d:1 o d 



fit 



(4) 



*t 



+ 4if D . D . Y (d*V$ 



d v Vt J )D^D 3 



where A d B = AdB — BdA; i,j are flavor indices; 
Vij = diag{(x>/\/2, uj/V2 : <ft} is the diagonal matrix con- 
taining id and (j) mesons (we omit the p and K* mesons); 
D* = {D*°,D*+,D*+) is the triplet, of vector D* 
mesons containing light antiquarks U, d and s, respec- 
tively. The chiral couplings g D , D , ^ , g D , D , v and f D , D , v 

are fixed as H-[I3|: 9 D * D , V = (3gv/V2, / D . D , V = 
m D ,Xg v /^/2, g D , D ,,^ = {m D *mj^)/{m D fj^), where 
fj, — 416 .4 MeV is the J/ip leptonic decay constant; 



gy ~ 5.8 and (3 w 0.9 are fixed using vector dominance; 
the parameter A = 0.56 GeV -1 is extracted by matching 
HHChPT to lattice QCD and light cone sum rules (see 
details in Q). The leading-order process relevant for the 
strong decays F(3940) -> J/ipu and F(4140) -> J/ip<j> is 
the diagram of Fig.l involving the vector mesons D* or 
D* in the loop. 

The coupling of the charged D* ± (D* ± ) mesons to pho- 
tons is generated by minimal substitution in the free La- 
grangian of these mesons. The corresponding electro- 
magnetic Lagrangian reads as: 

Com = eA a (g atJ D*-id^D*+ - D*~ id a D*+ + H.c 



+ e 2 D*-D*+ (A»A V - g^A a A c 



(5) 



This Lagrangian results in the two relevant diagrams dis- 
played in Figs. 2(a) and 2(b). In order to fulfill elec- 
tromagnetic gauge invariance the strong interaction La- 
grangian Cy also has to be modified. As outlined 
in Ref. and extensively used in Refs. 13, [l7j ]. 

each charged constituent meson field i/ ± in Cy is 
multiplied by the gauge field exponential: H ± (y) — > 
e *ieI(y,x,P)H±(y), where I(x,y,P) = $*dz^{z). Ex- 
panding e TieII - y ' x ' P ' > up to second order in the electromag- 
netic field, the two additional diagrams of Figs. 2(c) and 
2(d) are generated, which are necessary to guarantee full 
gauge invariance. The contribution of these additional 
processes is significantly suppressed (of the order of a few 
percent) compared to the leading diagram of Fig. 2(a). 

The invariant matrix elements of the strong and two- 
photon transitions (when all initial and final particles are 
on their mass shell) are given by: 



AV(F^ J/i>V) = 9„u9 YJ ^ 
M^(Y ->■ 77) = (g^q^ 



+ V 2fl v lu f Y J ^ 



(6) 



where t>i(<;i) and i>2(<72) are the 4- velocities (momenta) 
of J,/, and V. The effective strong couplings g YJ v an( i 
f YJ v have dimension of mass, while the electromagnetic 
coupling g Y has dimension of inverse mass. The matrix 
element of the two-photon transition has a full gauge- 
invariant structure. The constants q v , ,, and f v , ,, 
are products of the coupling gy, the chiral couplings in 
Eq. ([4} and the generic D* meson loop structure integral 
(see Fig.l). In terms of these effective couplings g YJ v , 
f YJ v and g Yjy the corresponding decay widths are cal- 
culated according to the expressions: 

3P* 



r(Y -> J/ipv) 

T(Y - 77) 



8nm Y 



2 

9yj,,,v 



(l + (3 + 2wrf3 + 3r 2 /3 2 ) 



2 S 2 



(7) 



where 
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FIG. 1: The diagram describes the Y — » J/tpV decay. 
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FIG. 2: The diagrams shows contributions to the Y — > 77 
decay. 



and a is the fine structure constant. Here P* is the cor- 
responding three-momentum of the decay products. 

Our numerical results for the quantities characterizing 
the strong J/ipV (V = oj,4>) and radiative two-gamma 
decays of F(3940) and F(4140) are contained in Table|U 
For the masses of the Y states we use the values extracted 
by the BABAR @ and the CDF [l[ Collaborations. The 
error bars correspond to the ones of the experimental 
mass values of the Y states. 

The predictions for the couplings <?y of the Y states 
to their meson constituents are consistent with a trivial 
estimate using the Weinberg formula. It was originally 
derived for the deuteron as based on the compositencss 



condition [13] with g^ 



3/4 1/4 



This formula 



TABLE I: Decay properties of F(3940) and Y~(4140) states 



Quantity 


Y(3940) 


Y(4140) 


gy, GeV 


14.08 ± 0.30 


13.20 ± 0.26 


SV^Vi GeV 


1.72 ± 0.03 


1.46 ± 0.03 


W>GeV 


1.64 ± 0.01 


1.84 ± 0.01 


r(y -> J/w), Mev 


5.47 ± 0.34 


3.26 ± 0.21 


g Y „ x 10 2 , GeV" 1 


1.15 ±0.01 


1.46 ±0.01 


T(Y -> 77), keV 


0.33 ± 0.01 


0.63 ± 0.01 


T(Y -» J/ipV) ° 


0.61 ±0.06 


1.93 ±0.16 
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represents the leading term of an expansion in powers 
of the binding energy e. Note that this expression can 
be obtained in the local limit (i.e. the vertex function 
approaches the limit $(y 2 ) — > <5 4 (j/)) and when the lon- 
gitudinal part k^k v /m 2 D , of the constituent vector me- 
son propagator is neglected. The numerical results for 

9Y(3940) = 9 ' 16 GeV and #Y(4140) = 8 ' 91 GeV &ie 1U S°° d 

agreement with nonlocal results of <?r(394o) = 14.08 GeV 
and gy(4i4o) = 13.20 GeV. 

The predictions of r(V(3940) -► J/ipu) = 5.47 MeV 
and r(V(4140) -> J/ip<j>) = 3.26 MeV for the observed 
decay modes are sizable and fully consistent with the 
upper limits set by present data on the total widths. 
The result for r(V(3940) — > J/ipco) is also consistent 
with the lower limit of about 1 MeV [4j. Values of 
a few MeV for these decay widths naturally arise in 
the hadronic molecule interpretation of the V(3940) and 
V(4140), whereas in a conventional charmonium inter- 
pretation the J/ipV decays are strongly suppressed by 
the Okubo, Zweig and Iizuka rule [4J. In addition to 
the possibility of binding the D*D* and D* + D*~ sys- 
tems [7[ , present results on the J/ ipV decays give further 
strong support to the interpretation of the Y states as 
heavy hadron molecules. Further tests of the presented 
scenario concern the two-photon decay widths, which we 
predict to be of the order of 1 keV. 

Finally we also test the J PC = 2 ++ assignment not yet 
ruled out experimentally. The coupling of the molecular 
tensor field Y^-ij to the two-meson constituent current 
Jy" is set up as 



C Yt {x) = gY T Y^ ;ij (x) I d 4 y<S>(y 2 ) J^{x,y) 



Proceeding as outlined before we obtain 



(8) 



r(V (3940) -> J/ipw) = 7.48 ± 0.27 MeV , 
r(V(4140) -> J/i/j4>) = 4.41 ± 0.16 MeV , 
r(V(3940) -> 77) = 0.27 ± 0.01 keV , 
r(V(4140) -> 77) = 0.50 ±0.01 keV . (9) 

Since the results for the strong J/ip decays are quite sim- 
ilar to the ++ case, a 2 ++ scenario cannot be ruled out 
and is also consistent within a molecular interpretation 
of the Y states. 

A full interpretation of the V(3940) and V(4140) states 
requires: i) an experimental determination of the J PC 
quantum numbers, ii) a consistent and hopefully con- 
verging study of binding mechanisms in the D* s ^D* s ^ 
systems and iii) theory and experiment to consider the 
open charm decay modes, such as DD, DD*, DD*j, etc., 
which are also naturally fed in a charmonium picture. Ul- 
timately, only a full understanding of the decay patterns 
of the V(3940) and V(4140) can lead to a unique struc- 
ture interpretation, yet present results clearly support 
the notion of the establishment of hadronic molecules in 
the meson spectrum. 



After submission of this manuscript, calculations both 
in the potential model approach [19j | and in QCD sum 
rules [20l | were presented, which support the original 
claim that the D*D* system binds for J PC = ++ , hence 
give further support to the interpretation presented here. 
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